Photonic band structure of highly deformable, self-assembling systems 
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We calculate the photonic band structure at normal incidence of highly deformable, self-assembling 
systems - cholesteric elastomers subjected to external stress. Cholesterics display brilliant reflec- 
tion and lasing owing to gaps in their photonic band structure. The band structure of cholesteric 
elastomers varies sensitively with strain, showing new gaps opening up and shifting in frequency. 
A novel prediction of a total band gap is made, and is expected to occur in the vicinity of the 
previously observed de Vries bandgap, which is only for one polarisation. 



PACS numbers: 61.30.-v, 81.40.Jj, 42.70.Qs 

Photonic band-gap (PBG) materials offer a new ap- 
proach to the manipulation of light that depends on the 
structure rather than the atomic or molecular properties 
of materials. These materials have two unique properties 
which has spurred interest in their design, namely the lo- 
calization of light JL| and modification of the spontaneous 
emission spectrum from atoms and molecules |2[] . Several 
approaches have been taken to manufacture PBG materi- 
als. Yablonovitch constructed an fee photonic crystal by 
drilling holes into a dielectric medium B. Later, Ozbay 
and co-workers designed a picket fence structure which is 
assembled by stacking two-dimensional layers 0. 

Recently, there has been an increased interest in self- 
assembling PBG systems due to their relative ease of 
manufacture for operation at optical and near-infrared 
wavelengths. Several examples include air holes in a tita- 
nia matrix jq| , copolymer-homopolymer films which form 
lamellar structures J6| , thin films of PMMA infilled with 
SnS 2 @, and cholesteric liquid crystals (CLC's) §Q. 

One of the most promising applications of pho- 
tonic band-gap materials is in low-threshold lasing. 
Yablonovitch first predicted that the lasing thresh- 
old would be decreased by introducing a defect into an 
otherwise perfect photonic material. Since spontaneous 
emission is suppressed in the bulk, excitation would not 
be drained by any emission into non-lasing modes. Such 
low-threshold lasing has recently been observed in two- 
dimensional photonic crystals |12[ . Alternatively, one can 
design lasers that take advantage of the enhanced dwell 
time associated with the band edge divergence of the 
density of states 13 1. Experimentally, this band-edge 
lasing has been observed in CLC's |1Q] and cholesteric 
elastomers (CE's) Q. 

A CLC has local orientational ordering along a direc- 
tor n, which rotates as a periodic function of distance 
along the pitch axis z. The director of an ideal CLC ad- 
vances uniformly, tracing out a helix of pitch pq. The 
pitch can be adjusted to match the wavelength of visible 
light, whereupon a number of spectacular optical effects 
are observed experimentally and explained theoretically 
P,P|. In particular, in experiments conducted at nor- 
mal incidence, circularly polarized light which twists in 
the same sense as the helix is reflected with its original 



polarization, while circularly polarized light that twists 
in the opposite sense is transmitted unchanged. Normal 
incidence has been of prime concern since the optical re- 
sponse of such twisting nematic media is the basis of 
liquid crystal (LC) display technology. A CLC can be 
considered locally uniaxial, with a dielectric permittiv- 
ity en along n and ej_ perpendicular to n. By solving 
Maxwell's equations in a rotating frame, de Vries found 
a single band gap in the photonic structure of an ideal 
CLC at normal incidence g . 

The calculations we present on CE's point to new phe- 
nomena and new applications, not possible in existing 
photonics and hitherto unsuspected in the liquid crystal 
field. For instance, we find multiple gaps, some not at 
the zone edges, in contrast to classical CLC's. We also 
observe gaps for light of the opposite handedness to the 
underlying helix, again unexpected in classical CLC sys- 
tems. At some points the gaps for both polarizations 
overlap, giving a total gap of significance when polari- 
sation control is required. Our systems are highly de- 
formable (to many 100s%) and we shall find shifts in 
the (developing) band structure that can be large. Exis- 
tent photonic media typically have piecewise variation of 
an isotropic refractive index in going between a matrix 
and its inclusions. By contrast, CLC's have a continuous 
variation of the principal axes of birefringence. Polari- 
sation effects are thus very subtle and become more so 
for oblique incidence, which we consider in greater detail 
elsewhere. Control of polarisation is at the heart of LC 
and optical devices; we thus view this work as a first step 
toward new classes of photonic solids with deformable, 
tunable band structures. 

CE's can be made by crosslinking cholesteric poly- 
mer liquid crystals |Q. Defect-free monodomain rubber 
strips tens of centimeters long display spectacular optical 
effects, viz. large changes in the frequencies of reflection 
and lasing mW in response to imposed mechanical strains 
that couple to director orientation (fig. |l|). These strips 
can be thick and are oriented not by surface anchoring 
as in liquids, but by interaction between local directors 
and the rubber matrix. 

Elongations A = X xx , applied perpendicularly to the 
pitch axis (see fig. B) are predicted to coarsen the initially 
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FIG. 1. A CE showing the initially helical director distri- 
bution, cj)o(z). Elongation \ xx is applied perpendicular to the 
pitch axis which contracts by \ zz . The dielectric tensor is 
represented in its local principal frame by ey and e±. 



helical director structure, given by tfio = q$z (where <j> 
is the angle the director makes with the x-axis), to one 
dominated by regions of slowly varying angles, separated 
by increasingly sharp twist walls |L5| : see fig. g. 

At a critical A = A c , the walls become thermodynam- 
ically unstable and the director experiences periodic os- 
cillations about 4> — which diminish with increasing A. 
There are attendant contractions perpendicular to the 
stretch. The pitch shrinks affmely pdj | with the matrix. 
In the small stretching limit (A — > 1), the pitch varies as 
P = PoA -2 / 7 and the first reciprocal lattice vector goes 
as q = <7oA 2 ' 7 . Thus, the band structure changes upon 
extension because of two factors: the dilation of the recip- 
rocal space and the change in the modulation character 
of the dielectric tensor along the pitch axis, e(z). These 
changes can be very large. 

Theory. — Maxwell's equations yield flqj : 

H = Vx [gizy 1 (V x H)] . (1) 

Consider normal incidence, along the z axis. The mag- 
netic field H is transverse and exists wholly in the xy 
plane. We thus suppress the z components in the inverse 



dielectric tensor, given by 





-i 



in its principal 



frame (oriented at angle cf>(z) to x, see fig. pi), and V is 
only (d/dz) z. 

We apply Bloch's theorem to decompose H into plane 
wave components mM, so that 



G,l 



(2) 



where the unit vectors are 7 = {1,2}, e% = x and 
&2 = y, and the reciprocal lattice vector G = 2nqz, 
for n integer. This procedure yields a matrix equation 
that reduces to a dimensionless form: lengths trans- 
form according to z — ► z = zq/2n, wave vectors go as 
k — > k = k/q, reciprocal lattice vectors G — > 2n, fre- 
quencies go as uj — > u> — uj/(cqy/a) and e _1 — > e _1 /a 
where a = §(— + — ). This reduction is important for 
the proper interpretation of the shifts of band structure 
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FIG. 2. Modification of director distribution, <f>(z) by ei- 
ther mechanical strain 8 (CE) or by DC magnetic fields h 
(LC). The dielectric tensor's principal frame follows 4>(z)\ dis- 
tortions to the helix then induce changes in the band struc- 
ture. 

with elongation in figs. & 0. Since q = q A 2 / 7 changes 
with A, so do u) and k. 
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then assumes the form A, , , .,h( n ~iV 

=(n7),(n7)' " n 



Eq. ® 

u} 2 h/ ni y A thus determines the photonic band struc- 
ture of a CLC. It depends on the reduced inverse dielec- 
tric tensor at arbitrary z and thus angle <f> = <f> (z): 



e -1 (z) = l-a (cos20)ct + (sin 2(f)) a 



(3) 



where a = (eii — ej.)/(eii + ej_) follows the notation of de 
Vries |j, and the a. are the Pauli spin matrices. 

One can then show that in 7 space, A is given by 



(4) 



A k , = (k + n) 2 S nn >l 

=n,ri y 1 n > n = 

+ a(k + n)(k + n') [c n >- n g_ z + s„<_ 

where s n /_„ and c n >^ n are the Fourier coeffi- 
cients of sin(20) and cos(20), respectively: s n = 
J dz sm[2(j)(z)] exp(— Airinz), and c n similarly. 

The undeformed director angles 0(5) are 0o = 27rz. On 
an re-strain A, accompanied by relaxation Aj,j,(A) assumed 
to be uniform and determined by energy-minimisation, 
the principal frame orientation is given by fig] 



tan 20 = 



2AAj,j,(r — 1) sin47rz 



(r - 1)(A 2 + A 2 ) cos4^z + (r + 1)(A 2 - A 2 ) ' 



where r is the shape anisotropy of the polymers under- 
lying the nematic phase. See fig. |2| for <ft(z) for various 
S = X — 1 . From tan 20, one can easily obtain sin 20 and 
cos 20, and thus, s n /_ n and <v_ n . 

Numerical diagonalisation of the matrix A, ,. . ,., at 

° =(nA),(nA)' 

a range of k yields a dispersion relation ui(k), along with 
eigenvectors giving the character of each solution. In 
general, the eigenvectors are elliptically polarised inside 
the CLC medium, with semi-major and semi-minor axes 
corresponding at each point to the local principal axes 
of the dielectric tensor, and nearly circularly polarised in 
vacuo S. We take ey = 3 and a depressed value e± = 1.2 
throughout, simply for readability. 



At small k <C 1, the dispersion relation for an ideal 
CLC (i.e., A = X yy — 1) is linear, corresponding to non- 
dispersive waves, with a simple effective refractive index 
m = \/(e± + e ||)/2j suggesting that both modes effec- 
tively experience the same, homogeneous medium at long 
wavelengths. This small k behaviour is initially retained 
in the strain-modified band structures, see fig. 0. 

At k = 1, the branch whose polarisation rotates in the 
same sense as the helix in vacuo develops the de Vries gap 
M . The eigenmodes of this branch at the zone bound- 
ary are linearly polarised inside the CLC medium. The 
lower band's electric vector points along n(r), the up- 
per band's perpendicularly to n(r), in the xy plane. The 
other branch, however, cannot split analogously, since its 
polarisation in vacuo rotates in a sense opposite to the 
helix. This is qualitatively like the major gap in the dis- 
torted band structure, marked with dots at k = 1 in fig. [|. 
Furthermore, no gaps are observed for k > 1, because 
band gaps are created only when degenerate energy states 
are linked by non-zero matrix elements. Since in the de 
Vries case, there are only two harmonic components of 
equal and opposite frequencies, only the matrix elements 
linking the two lowest energy states on either side of the 
first Brillouin zone boundary are non- vanishing. 

We now stretch CE's with r — 1.9 for definiteness, 
which gives a A c « r 2 / 7 rs 1.2 |[J||. Fig. || shows the 
dispersion relation for an elongation A = 1.1 < A c . Since 
the first Brillouin zone boundary is at k = 1, band gaps 
may occur at k — nqo\@ , for integer n ^ 0, with (3 = 2/7. 
This corresponds to a shift in colour and lasing frequency 
[fill toward the ultraviolet. 
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FIG. 3. For A = 1.1 < A c , two differences from the de 
Vries case are observed: (1) additional bandgaps are created 
at higher zone boundaries for the RP branch and (2) bandgaps 
are also observed for the previously uninteresting LP branch, 
albeit much smaller than the RP gaps. The single gap in the 
de Vries case is approximately that marked by dots. 
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FIG. 4. For A = 1.3 > A c , substantial divergence from 
the de Vries dispersion relation is observed. A full band gap 
away from the Brillouin zone boundary is observed, as well as 
several anti-crossings between branches. 

For A > A c , there is a qualitative change in the be- 
haviour of the director, 4>(z) (see fig. g), and thus a qual- 
itative change in the band structure. Additionally, the 
scaling behaviour of X ZZ (X) = X~^ changes ]15| from 
j3 = 2/7, in the limit of small stretching (A ~ 1), to 
(3 = 1/2 for A > A c , the classical exponent for isotropic 
CE's. Fig. [| shows the dispersion relation for a stretch 
A = 1.3 > A c . 

We now analyse the gap structures that open up in the 
stretched case, A > 1. The elastic strain, S = A — 1, is 
the perturbation parameter modifying the perfect helical 
structure. Whereas before, c±\ — 1/2, and c„ = other- 
wise, we now have non-zero values for c± n which scale as 
5 n - 1 for n > 1 and S < 1; c ±1 = 1 - O (d 2 ) and c ~ 6. 
Applying degenerate perturbation theory, we eliminate 
all matrix elements in A™° , , ., except for those linking 

the degenerate energy states, and predict that the gaps 
for the interesting polarisation in the de Vries case will 
scale with the magnitude of the off-diagonal elements, 
given by c na . The 1 order (de Vries) gap will be ap- 
proximately constant, the 2 order gap (first new gap) 
will grow as S, the 3 order gap will grow as S 2 , and 
so on: see fig. pi More precisely, at the I s order gap, 

lu 2 = 1 ± a (ci ± \/c\ — s\ 1 , since s = 0. As a result, 
we also predict that the opposite polarisation will now 
have a non-zero gap that scales as 5 2 . That tells us that 
a full photonic band gap will be created at k — 1 as we 
perturb our helix by transverse elongation, and that its 
size will scale as S 2 , the size of the much smaller gap of 
opposite polarisation. See figs. |^ and ||. 

Of experimental interest is the in vacuo wavelength, 
A, of the light corresponding to a given Co on the disper- 
sion relation, particularly at the gaps. The definitions 
below (||) give A = p /(Lby/aX 2 ^ 7 ). Pitches po typically 
give a band in the visible so the initial wavelengths are 
Ao = Po/\/a ~ 500nm at Q = A = 1, which allows us 
to write A = Ao/((I»A 2 / 7 ). Likewise the first order de 
Vries gap is given by AA w Ao/A 2//7 a. The higher or- 



der gaps of the same polarization will have widths of 
AA„ sa CAoS n ~ 1 /(nA 2 / 7 ), where C is a pre-factor of or- 
der unity that will depend on r and a. For example, the 
second order gap in a rubber with r = 1.9, a = 0.43, 
and 8 = 0.1, the gap will be AA 2 sw 0.045A . For 
Ao ~ 800nm, that implies a stop band for the light with 
a circular polarization of the same sense as the helix will 
be observed for A = 362 nm to A = 398 nm. 

Finally, a note about oblique incidence in CE's. By 
symmetry, we expect that the magnetic field must have 
the same magnitude at all points for a given z, and only 
differ by a phase. That lets us generalize the H vector 
for normal incidence, eq. (J2|), by k — > k' = k±p + k\\Z 



and e 7 — > e, G y Our preliminary findings indicate that 
for a constant k' , the stop bands shift upwards as we 
increase the angle of incidence from zero, which implies 
that refraction out of a normally incident beam path is 
forbidden for modes just above the stop band. That ob- 
servation provides a mechanism to explain the spatial 
coherence of the light produced by dye-doped pumped 
lasers based on CLC's and CE's [0Q. 
Cholesteric Liquids. — We apply an external magnetic 
field H along the y direction. A CLC has an anisotropic 
susceptibility Xa = X\\~ X±- The helix untwists (increas- 
ing the period) and coarsens [|l7J until the energy gain 
from aligning with the field balances the Frank penalty 
for deviations from the original structure. The coars- 
ening of the director orientation is illustrated in fig. g, 
with the z-coordinate being reduced to z by the length- 
ening period. At a critical field H c = (irqo/2)^K22/Xa, 
where K22 is the Frank twist elastic constant, the period 
diverges logarithmically as the entire sample aligns with 
the external field |E7j . For typical cholesteric liquids with 
a pitch of 20 fjxa., H c = 15, 000 G and E c = 50 sV/cm 
Pj . Differences from the case of a CE under strain are 
detailed in |Q. 

The optical implications of coarsening were investi- 
gated by Meyer p8[ . The Fourier coefficients describing 
the twist of the cholesteric liquid are found to scale just 
like the coefficients describing CE's for A < A c . One can 
use the results for the CE case for a cholesteric liquid un- 
der the transformation 8 — > h 2 . The dispersion relations 
are qualitatively the same as in fig. 0. The width and 
scaling of the gaps created resemble those of fig. M, but 
with all effects ending abruptly at h = 1. 

Conclusions. — An entirely new type of photonic mate- 
rial has been described and characterized. Not only is it 
self-assembling and easily available as large, defect-free 
single crystals, but it is highly deformable. Earlier de- 
scriptions |TJj,[l7| °f hs modified periodic dielectric struc- 
ture have been used as the basis for calculating its band 
structure. New gaps arise and their widths scale in a 
well-understood fashion with the stretch applied to the 
material or the strength of the external field. The midgap 
frequencies shift position by large amounts comparable 
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FIG. 5. Scaling of the first to fourth gap sizes of a CE with 
8. The points represent numerical data, the straight lines, 
predictions from perturbation theory (assuming a scaling form 
for Cn and s n ). 

to their initial values. 
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